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In the present work we show that planetary mean distances can be calculated with the help of 
a Schrodinger-type diffusion equation. The obtained results are shown to agree with the observed 
orbits of all the planets and of the asteroid belt in the solar system, with only three empty states. 
Furthermore, the equation solutions predict a fundamental orbit at 0.05 AU from solar-type stars, a 
result confirmed by recent discoveries. In contrast to other similar approaches previously presented 
in the literature, we take into account the flatness of the solar system, by considering the flat 
solutions of the Schrodinger-type equation. The model has just one input parameter, given by the 
mean distance of Mercury. 



I. INTRODUCTION 

In recent years, some authors Q have suggested a 
quantum-like approach to calculate the planetary orbits 
in our solar system, which has led to several impressive 
results. Besides obtaining the observed orbits of all the 
planets and of the asteroid belt, such models have made 
predictions subsequently confirmed by observations. 

For example, the existence of asteroids orbiting be- 
tween Uranus and Neptune, with orbit radii around 24.77 
AU p] , was confirmed by recent discoveries . Further- 
more, the prediction of two intramercurial orbits with 
radii around 0.05 AU and 0.18 AU 0,0], absents of our 
solar system, has been verified in several extra-solar plan- 
etary systems during the last years [T^ . 

In spite of these and others positive results, these ap- 
proaches seem to be very speculative, for the quantization 
of macroscopic systems is something outside the scope of 
our known physics. Nevertheless, some possible origins 
for such effects have been outlined, some of them on the 
basis of quite orthodox theories. 

A large scale quantization can be based, for instance, 
on an extension of the ordinary commutation rules, in 
order to recover the equivalence principle in the con- 
text of quantum mechanics [3, Il3j . The general com- 
mutation rules derived in this way predict two scales of 
quantization: the usual microscopic one, and the scale in 
which the gravitational interactions are dominant. The 
presence of quantization at so different scales can be re- 
lated as well to self-similarity concepts 0. And it is 
also characteristic of fractal space-time approaches 0, 
in which the particle trajectories are non-diffcrcntiable 
at very short and very large scales. 

The concept of non-differentiability of space-time at 
short scales has been used by some authors in trying to 
derive the Schrodinger equation in the context of a classi- 
cal framework. For example, in the stochastic mechanics 
of Nelson [l4|, the Schrodinger equation is obtained as 
a classical diffusion equation, with help of the hypothe- 
sis that any particle in the empty space, under the influ- 
ence of any interaction field, is also subject to a universal 



Brownian motion without viscosity. Under this hypoth- 
esis, a Schrodinger-type equation follows from classical 
dynamics, with the Planck constant being related to a 
diffusion coefficient. 

The main problem of this kind of derivation of the 
Schrodinger equation (as well as of other descriptions 
based on diffusion processes or fluid dynamics [l5j|l is 
looking for a convincing physical origin for that universal 
Brownian motion, although some possibilities have been 
outlined, based for instance on quantum fluctuations on 
cosmic scale |16| , or on the quantum nature of space-time 
in quantum gravity theories |17j . 

Nevertheless, the important point in Nelson's work is 
that a diffusion process can be described in terms of a 
Schrodinger-type equation. In this context, the possibil- 
ity of describing a classical process like the formation of 
solar system in terms of quantum mechanics can be se- 
riously considered. The use of such an approach is also 
suggested by the chaotic behaviour of the solar system 
during its formation and evolution |l8j-|2lj. which im- 
plies the non-differentiability of trajectories for large time 
scales 0. 

Despite all these conjectures, the physical principles 
behind the quantization of large structures cannot be 
considered completely understood. In the present work, 
our aim is not to solve this problem, but just to explore a 
little more the results we can obtain from a quantum-like 
approach. Our goal is to introduce a new ingredient not 
considered up to now: the flatness of the solar system. 
For this purpose, we will look for flat solutions of the 
Schrodinger-type equation, in contrast to the spherical 
solutions usually considered in the previous papers. 



II. PRELIMINARY REMARKS 

Applying a Bohr-like "quantization" rule to the orbital 
angular momentum of a planet in circular orbit around 
the Sun, L = mvr = ng*/2ir, and using the Newtonian 
law for the orbital velocity, v 2 — GM/r, one of the au- 
thors has derived the following expression for the radii 
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of the orbits: 
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where n is an integer, G is the gravitational constant, M 
is the solar mass, m is the planet mass and g* stands 
for a quantum of action, playing the role of a re-scaled 
Planck constant. 

In the case n — 1, equalizing Q to the observed value 
of the Mercury orbital radius, using the observed values 
for G and M, and taking m = 2.1 x 10 26 kg (the av- 
erage mass of the planets of the solar system), we find 
g* = 3.6 x 10 42 J.s. This value coincides to a re-scaled 
quantum of action derived on the basis of self-similarity 
considerations 0. 

Now, equation Q can be written as r — n 2 ri, where 
r*i is the Mercury orbital radius. Taking for n the next 
integer values, we can obtain a sequence of values that 
assures very well the observed values of orbital radii in 
our solar system pj. For n — 2 we obtain the orbit of 
Mars; for n — 3 the orbit of Camilla (an asteroid located 
in the outer region of the asteroid belt between Mars and 
Jupiter); the values n = 4 and n — 5 give the orbits of 
Jupiter and Saturn, respectively; for n = 6 we have the 
orbital radius of Chiron (an asteroid between Saturn and 
Uranus); and, for n — 7, 9 and 10, the orbits of Uranus, 
Neptune and Pluto follow, in this order. The deviations 
from the observed values are 19% for Jupiter, 4% for 
Neptune and less than 2% in the other cases. 

The value n = 8 predicts an orbit between Uranus 
and Neptune, with radius 24.8 AU. In reference pj, sub- 
mitted for publication in October 1994, this orbit was 
considered empty. Surprisingly enough, the value above 
is in very good agreement with subsequent discoveries 
of two asteroids located precisely at this radial position, 
named 1993 HA 2 and 1995 DW 2 [HI These discov- 
eries were published in 1995, that is, one year after the 
paper submission, showing the predict power of this kind 
of model. 

The most important drawback of the model is, as al- 
ready noted by the reader, the absence of Venus and 
Earth in the above sequence. To include these planets, 
and also the asteroid Vesta, located in the inner region 
of the asteroid belt, it has been taken into account an ad 
hoc second quantum number n' , running from to n 0, 
such that relation is generalized to 
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For n = n' we recover the original relation and results. 

Another possibility was provided by Agnese and Festa 
0, taking for n = 1, instead of the Mercury radius, the 
value r\ = 0.04 AU, which corresponds to g* = 1.2 x 
10 42 J.s. Now, for n running from 3 to 6, we obtain 
the orbital radius of Mercury, Venus, Earth and Mars; 
for n — 8 we have the asteroid Ceres; and, for n = 11, 
15, 21, 26 and 30, the orbital radius of Jupiter, Saturn, 



Uranus, Neptune and Pluto follow, in agreement with 
the observed values. The problem here, as we can see, 
is that we fall into a lot of empty orbital positions, that 
is, radius values predicted by the formula r = n 2 ri, but 
not occupied by any observed body, particularly for large 
values of n. Nevertheless, the prediction of a fundamental 
radius at 0.04 AU would be shown an important result. 

Another interesting contribution by Agnese and Festa 
is the introduction of the re-scaled fine structure constant 
ct g = 2TrGMm/g*c, where c is the light velocity. In this 
way one can rewrite equation Q as 
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Using for m the average mass of planets of our solar sys- 
tem and for g* the value obtained above, we obtain for 
a g the value used by Agnese and Festa, a~ 1 = 2.1 x 10 3 . 

This quantum approach also works for other, recently 
discovered, planetary systems, as shown by Nottale 0,0 
and also by Agnese and Festa 0, @- In Nottale's ap- 
proach the periods of planets revolution are quantized 
as 



T = 27rGMn 3 /wo, 



(4) 



where ujq = 144 km/s is a characteristic velocity. Using 
T = 2irr/v, v 2 = GM/r and r given by ®, we derive 
= ctgC- Using ujq — 144 km/s, we obtain a^ 1 — 
2.1 x 10 3 , which again is in accordance with the value 
obtained by Agnesa and Festa. 

Therefore, Nottale also predicts a fundamental radius 
given by r\ ~ 0.04 AU. Several extra solar planets re- 
cently discovered lie at this distance from their star ^1 • 
This agreement between a fundamental radius derived 
years ago and subsequent observations shows once again 
the predict power of such models. 

Some of the above results come from a semi-classical 
treatment of the problem, based essentially on the appli- 
cation of a Bohr-like quantization condition. In Nottale's 
works we can find a more rigorous approach, based on a 
Schrodinger-type wave equation. In this way he can ac- 
commodate all the planets and the asteroid belt of our 
solar system, with only a few empty orbits |3J| . To achieve 
this goal he treats the inner and outer systems separately, 
with two distinct values of ujq. 

In the next section, we propose to show that flat so- 
lutions of the Schrodinger-type equation can furnish im- 
portant information on the mean planetary radii. Our 
aim is to enrich the previous studies, including a new 
ingredient not considered up to now: the flatness of the 
solar system and of the original disk from which it was 
originated. 



III. THE FLAT SOLUTIONS 

In this section we present the mathematical results 
obtained from the solution of a Schrodinger-type equa- 
tion involving an attractive central field. As already 
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mentioned, the orbits of planets and asteroids will be 
considered approximately in the same plane, and the 3- 
dimensional equation will be solved taking z = 0. 

Let us consider a body of mass m moving around an- 
other body of mass M , under an attractive central field 
with potential V(r), where r is the axial distance in po- 
lar coordinates. The Schrodinger-type equation is given 
then by 



d 2 ip ldip 1 d 2 ifj 



2p \dr 2 



di 



r 2 QQ2 



+ V(r)tp = E%1). (5) 



The term in parenthesis is A-0, where A is the Laplace 
operator in polar coordinates (r, 0) , which is applied to 
the wavefunction tp of the body of mass m. The param- 
eter E stands for the total energy of the system; p is the 
reduced mass mM/(m + M); and g* is a constant, as 
previously described. 

As long as the potential V is a function of the radial 
variable only, we may look for a solution using separation 
of variables, 

lKr,0) = /(r)0(0). (6) 
Replacing © in JSJ) and dividing by /(r)6(0), we get 
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2p f(r) 



d 2 f(r) | ldf(r) 
dr 2 r dr 



g* 2 1 d 2 0{8) 



+ E-V{r) = 0, 



2pr 2 6(0) d9 2 
or, equivalently, 

1 d 2 Q(9) r 2 \d 2 f{r) ldf(r) 
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Since the term on the left hand side depends only on 
and the term on the right hand side depends only on r, 
both terms must be equal to a constant which we denote 
by — £ 2 . Therefore, we obtain two ordinary differential 
equations: 



6"(0) - -£ 2 G{6), 



(9) 



and 



f"(r) + l -f'(r) + {-*+?P[E- V(r)} \ f(r) = 0. 



Equation has a solution 

6(0) = e m . 
If we assume the boundary condition 

8(0) - 0(27r), 
we end up with \l\ = 0, 1, 2, 3 . . .; I is an integer. 



(10) 
(11) 

(12) 



Now we apply a rescaling in (|10|) by 

p = 2/3r, P>0, 
where 1 = -2pE / g* 2 . Further more we define 
_ fxGMm 



(13) 



(14) 



and, by using V(r) — —GMm/r, equation (|TU|l becomes 



}'(p) + -Hp) + 
p 
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+ - f( P ) = 0. (15) 



If we consider f(p) = u(p) in 11511 we obtain 

Vp 



n (I 2 -1/4) 
P P 2 



i(p) = 0. (16) 



For this equation, p = is a regular singular point and 
p = oo is an irregular singular point. Equation l|16|) is a 
confluent hypergeometric equation p2l | , referred as Whit- 
taker 's equation [23j. This equation has a regular solu- 
tion given by a hypergeometric series which converges if, 
and only if, 



±-±k, 
2 



k = 0, 1,2,3. 



(17) 



By definition, n > 0. Assuming that the solution u(p) 
satisfies the boundary condition 



lim u(p) = 0, 



and as u(0) must be finite, we have 



n = \t\ + - + k, 



fc = 0,l,2,3... 



So far we have the conditions 
|^| = 0,1,2,3... 
n = + + | \£\ + \- 

They can be regrouped in the form 

13 5 7 

U ~ 2'2'2'2"' 



(18) 



(19) 



(20) 



0,1,2,. 



1 

,n--. 



(21) 



For each pair n£, the solutions u n £ were obtained us- 
ing Maple V Release 4 Software (Waterloo Maple Inc.). 
As long as we obtain u n i(p), we may find the solution 
fntir) of the original equation 1|10|) . Once we know the 
collection of functions f n i(f) and Q^(r), their products 
generate the normalized solutions Vw of the Schrodinger- 
type equation. 
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Now we define 



Pnt(r)dr = / ip^nirdrdd, 
Jo 



which is equivalent to 



P n e(r)dr = r[f(r)} 2 dr. 



(22) 



(23) 



As in the case of atomic wavefunctions, the mean radius 
is obtained by 

/>oo />oo /•2-k 

r n t = / rP ne (r)dr= / / (ip^rip ne ) rdrd6 , (24) 
Jo Jo Jo 

which is equivalent to 



r n i = / [rf{r)fdr. 
Jo 



(25) 



In this context, we now return to the rescaling factor 
j3 and consider the mean distance of Mercury as r| . (If 
we associate Mercury with the first state ri , all results 
are inconsistent with the observed mean planetary dis- 
tances. The orbits of Venus, Earth, Mars and Jupiter, 
for example, definitely cannot be derived.) For n = 3/2 
and I = 0, we obtain for Eq. I)16|l the solution 



—u{p) = ci ( -e 
VP 



1-2P -L 0—5P 



e , 



(26) 



and so /(r) = c ie -^ r {2(3r - 1). Since / r[f(r)] 2 dr = 1 

Jo 



we have 



n = \/ g/3- 



From Eq. (|25|) we then get 



r |0 = 1 ^= 0.387 AU, 



(27) 



(28) 



where we have used the observed mean distance of Mer- 
cury. So, we obtain 



1 7 



5 0.387 3 



9.04/(3/2) AU- 1 . (29) 



Therefore, since (3 is proportional to 1/n (see Eq. H14H ). 
for each value of n we will consider 



0n = 9.04 (1/n) AU- 1 . 



(30) 



Its important to remark that the orbital radius of Mer- 
cury will be the only input parameter used in this model. 



IV. RESULTS AND DISCUSSIONS 

We denote by the triple (n, £, r) the resulting numbers 
n, £, and the corresponding mean planetary radius r n i (in 



AU) , calculated from equation (|25|l . The first triple found 
(1/2,0,0.055) does not correspond to any observed orbit 
in our solar system. As previously mentioned, Nottale 
jj| also found this radius as the fundamental level in his 
quantization law, and it was also obtained by Agnese and 
Festa as the first state in their quantization condition • 
As remarked in these two papers, the first extra-solar 
planet discovered around the solar- type 51 Pegasus lies 
precisely at 0.05 AU from its star. Other subsequently 
discovered companions of solar- type stars, which also fall 
around 0.05 AU, confirm this result [l^. 

The following two states (3/2,0,0.387) and 
(3/2,1,0.332) are associated with the orbit of Mer- 
cury. The triples (5/2,0,1.05) and (5/2,1,0.995) have 
very similar values for r which are equal to the mean 
planetary distance of Earth. The triple corresponding 
to Venus is (5/2,2,0.83). 

The set of triples (7/2,0,2.04), (7/2,1,1.99) and 
(7/2,2,1.82) does not correspond to any observed planet 
in this region of the solar system, but could be associ- 
ated with the asteroid Hungarias, at 1.94 AU. This radius 
value was also predicted by Nottale in his quantization 
law 0] and nearly corresponds to the radius found for a 
companion of the solar- type star 47 UMa B, lying at 2.12 
AU from this star (see 0, an d references therein) . 

The state (7/2,3,1.54) clearly stands for the orbit of 
Mars. The results obtained for n = 9/2 are (9/2,0,3.37), 
(9/2,1,3.32), (9/2,2,3.15), (9/2,3,2.88) and (9/2,4,2.49). 
They are in a very good correspondence with the asteroid 
belt in the solar system. The lower values 2.88 AU and 
2.49 AU correspond to the mean distances observed for 
the Ceres group (2.64 AU), the central peak of the belt, 
and for Vesta (2.36 AU), which delimits its interior ring. 
The upper values 3.32 AU and 3.37 AU correspond to 
the mean distance observed for Camilla (3.48 AU), which 
delimits the exterior ring. The intermediate value 3.15 
AU has a precise correspondence with the Hygeia group 
(3.16 AU), the main peak of the asteroid belt. 

For larger values of n, a great number of states is found. 
Nevertheless (as pointed out by Nottale et al. to ex- 
plain the circularity of the orbits in the solar system), 
a great number of states occupying the same region of 
space, most of them with large eccentricities, leads to a 
strong chaos and to the crossing of orbits, which on large 
time scales would generate the condensation of states on 
the observed, approximately circular orbits. Therefore, 
for n > 9/2 we will consider only the states with rota- 
tional symmetry, that is, with / = (see Eq. Ijlljl'). 

In this way, the states (11/2,0,5.03) and (15/2,0,9.34) 
can be associated with the orbits of Jupiter and Saturn, 
respectively. The states (17/2,0,12.0) and (19/2,0,15.0) 
give mean values around the orbit of the asteroid Chiron, 
distant 13.7 AU from the Sun. The state (21/2,0,18.3) 
can be associated with the orbit of Uranus, while the 
state (25/2,0,25.9) has a mean distance close to the or- 
bital radii of the already referred asteroids 1993 HA 2 
and 1995 DW 2. The next state, (27/2,0,30.2), clearly 
stands for the orbit of Neptune and, finally, the state 
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(31/2,0,39.9) can easily be associated with the orbit of 
Pluto. 

As one can see, all the planets in the solar system (as 
well as the asteroid belt) are fitted by the model. Apart 
the fundamental radius 0.055 AU (which should be ruled 
out for thermodynamical reasons @ ) , we have only three 
empty orbits, corresponding to the states (13/2,0,7.02), 
(23/2,0,21.9) and (29/2,0,34.9). In what concerns the 
occupied states, the deviations from the observed orbital 
radii are 15% for Venus, 10% for Chiron, 6% for Ceres and 
Vesta, 5% for Uranus and the asteroids HA2 and DW2, 
3% for Jupiter, Camilla and Hungarias, 2% for Saturn, 
1% for Mars and Pluto, and 0.5% for Earth, Neptune and 
Hygeia. 

Beyond the orbit of Pluto, we still have the states 
(33/2,0,45.2) and (35/2,0,50.8), which could correspond 
to other mass peaks in the Kuiper belt, a group of about 
60 trans- Nep tunian objects orbiting up to 50 AU from 
the Sun |2J, |25j . [After the completion of this work, it 
was announced the discovery of Quaoar, the biggest as- 
teroid discovered up to now in the Kuiper belt j2g. Its 
stable orbit has mean radius around 42 AU, which may 
correspond either to the state (31/2,0,39.9), associated to 
Pluto's peak, or to the state (33/2,0,45.2). The relative 
deviations are 5% and 8%, respectively] 

V. CONCLUSIONS 

The set of obtained results in calculating planetary or- 
bits in the present work, as well as those obtained by 
other authors, are enough to encourage further studies. 
In all these theoretical approaches, the existence of a 
fundamental radius around 0.05 AU is predicted, and 
several planets have been recently discovered orbiting at 
this radius in extra-solar systems. Furthermore, it was 
predicted the existence of asteroids between Uranus and 
Neptune pj, a prediction confirmed by observations. 

The Scrodinger-type model presented in this work 
seems to describe the observed mean distances from Mer- 
cury until the asteroid belt. For the outer system we have 
considered only states with rotational symmetry, obtain- 
ing good results as well, with only three empty states. All 
the results were obtained with just one input parameter, 
namely the orbital radius of Mercury. 

The reader may be asking how a set of planets and 
asteroids with so different masses can be described by just 
one equation. In other words, why do we need only one 
input parameter (the value of for Mercury) , if equation 
© depends explicitly on the reduced mass /i? 

The answer may reside on the fact that our 



Schrodinger-type equation is not a genuine quantum 
equation, but describes a kind of diffusion process. In this 
context, the re-scaled "Planck constant" present in (0 is 
related to the diffusion coefficient V by g* = 2mV |5lll4|. 
where to is the mass under diffusion. Furthermore, for 
any planet or asteroid gravitating around the Sun, the 
reduced mass coincides with its mass m. Therefore, 
remembering that V(r) and E are also proportional to 
to, it is easy to verify that JSJl, in fact, does not depend 
on m. 

Another way of verifying this result is to look at equa- 
tion 1141) . Using fi — to and g* = 2mT>, we obtain 
j3 = GAf/(42? 2 ?i). So, for a given diffusion coefficient 
and a given stellar mass, the parameter f3 depends just 
on n, as assumed in Eq. (|30|l . 

As a last remark, let us point out an observable differ- 
ence between the flat model presented here and spherical 
models. In the approaches by Nottale et al @ and by 
Agnese and Festa 0, there is a predicted orbit at 0.18 
AU from the Sun, corresponding to n = 2 (remember 
that, in their models, Mercury corresponds to n = 3). 
This orbit is also predicted through the Bohr-like model 
presented in , corresponding to the quantum numbers 
n = 1, n' = (see Eq. (J3|). As a little planet orbiting at 
the fundamental radius 0.05 AU is ruled out by thermo- 
dynamical constraints J5j, a body at 0.18 AU would be 
the only new body predicted to orbit near the Sun. 

As we have seen in the previous section, such an or- 
bit does not exist in the present model. It is excluded 
by the flatness of our solutions. Of course, this does 
not mean that it could not be observed in extra-solar 
systems. Actually, orbits with radius around this value 
were already observed This fact could suggest that, 
in contrast to our solar system, the corresponding plane- 
tary systems are not flat, which could indicate that they 
were originated in a different way. 
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